Evidence for forward scattering and coupling to acoustic phonon modes in high-Tc 

cuprate superconductors 



S. Johnstoni'4, I. M. Vishiki'2'3, W. S. Lce^'^ F. Schmitti'^, S. Uchida^, K. 
Fujita^, S. Ishida^, N. Nagaosa^'^, Z. X. and T. P. Devereaux^-^ 

^Stanford Institute for Materials and Energy Science, 
SLAG National Accelerator Laboratory and Stanford University, Stanford, CA 94305, USA 
^Geballe Laboratory for Advanced Materials, Stanford University, Stanford, CA 94305, USA 
^Department of Physics and Applied Physics, Stanford University, Stanford, CA 94305, USA 
*IFW Dresden, P.O. Box 27 01 16, D-01171 Dresden, Germany 
"Department of Physics, Graduate School of Science, 
University of Tokyo, Bunkyo-Ku, Tokyo 113-0033, Japan 
^Laboratory for Atomic and Solid State Physics, Department of Physics, 
Cornell University, Ithaca, New Yori: 14S53, USA 
''Department of Applied Physics, University of Tokyo, Bunkyo-ka, Tokyo 113-8656, Japan and 
^Cross-Correlated Materials Research Group (CMRG) and Correlated 
Electron Research Group (CERG), RIKEN-ASI, Wako 351-0198, Japan 
(Dated: January 10, 2011) 

Recent laser angle-resolved photocmission spectroscopy studies have established the presence of 
a new kink in the low-energy nodal dispersion of Bi2Sr2CaCu208+5 (Bi-2212). The energy scale 
(~ 8 — 15 meV) of this kink appears below the maximum of the superconducting gap Aq. Therefore 
it is difficult to interpret this feature in terms of the usual coupling to a sharp dispcrsionlcss mode. 
In this paper we examine electron-phonon coupling to the in-plane acoustic phonon branch arising 
from the modulation of the screened Coulomb potential. We demonstrate that such a coupling 
has a strong forward scattering peak, and as a consequence, a kink occurs in the dispersion at an 
energy scale shifted by the local gap A(k). In addition, considerations for the reduction of screening 
with underdoping naturally explains the observed doping dependence of the low-energy kink. These 
results point to a strong coupling to the acoustic branch which is peaked in the forward scattering 
direction and has important implications for transport and pairing in the high-Tc cuprates. 



PACS numbers: 71.38.-k,74.72.-h,74.25.Jb 

Understanding what controls the superconducting 
transition temperature T^. in the cuprates remains one 
of the major questions in condensed matter physics. The 
large variation in Tg across cuprate famihes remains un- 
explained despite the near similarity of the Cu02 plane 
and accompanying strong correlation physics that well 
describes the parent phases of these systems. Strong 
electron correlations are of course important in all the 
transition metal oxides, h\ii the cuprates arc special in 
the sense that they are the most metallic system among 
them. As a result, the cuprates are at the verge of lo- 
calization/delocalization and/or ionicity/covalency. This 
subtle balance or dichotomy may be key to the origin of 
the robust superconductivity, which requires both strong 
coupling and quantum coherence. 

This battleground has recently been highlighted in 
angle-resolved photocmission (ARPES), where polaronic 
lattice effects demonstrably alter the spectral function 
from a Lorentzian to Gaussian lineshape^'^ in undoped 
systems. The presence of dispersion "kinks" in doped 
cuprates have been (controversially) interpreted as cou- 
pling to lattice phonons, with stronger couplings than 
implied by density functional calciilations.^^'"' Finally, 
with the advancements of laser- based ARPES, a new 
energy scale has been observed in the nodal region of 
Bi2Sr2CaCu208+5 (Bi-2212) which manifests as a subtle 
kink in the band dispersion at an energy 8-15 meV.^^^° 

In this paper, we focus on the latter kink in under- 
doped Bi-2212 (Tc = 55K, UD55, p ~ 0.088) and present 



data and analysis which suggests that the interplay of 
Coulomb interactions and lattice effects as a consequence 
of ionicity/covalency dichotomy plays an important role 
in the physics of a doped Mott insulator. The experimen- 
tal setup is identical to Ref. 9. The energy resolution was 
3 mcV and no dcconvolution methods have been applied 
to the data.'' Whereas previous studies focused on the 
low-energy kink at the node, we have done experiments 
away from the node where the gap is nonzero. Since the 
opening of a gap will affect the momentum distribution 
curve (MDC) dispersion close to the gap energy, we chose 
a doping where the low-energy kink is very strong (Ref. 
9) so the two features can be easily distinguished. How- 
ever, we emphasize that the results presented here hold 
for other dopings in the undcrdopcd regime. 

The main experimental observations for the low-energy 
kink along the nodal direction are summarized in Fig. 1. 
As evident in the raw data (Fig. la), the kink mani- 
fests as a bend in the dispersion and an accompanying 
change in spectral intensity at an energy ^ 8-15 mcV in 
UD55 Bi-2212. This behavior is also evident in the MDC 
derived dispersion, which deviates from the linear disper- 
sion at the same energy scale (Figs. lb). For UD55, the 
slope of the dispersion at is considerably smaller than 
the slope of the dispersion between 30-40 meV (umid) and 
the low-energy kink is unmistakable. This observation 
has already been discussed in the context of the univer- 
sal nodal Fermi velocity,® suggesting the importance of 
this new energy scale in the low-energy nodal physics 
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FIG. 1: (Color online) Nodal low-energy kink in Bi-2212 
UD55. (a) Color image plot of raw data. The 70 meV and 
low-energy kinks are marked by arrows and both can be seen 
by eye in the raw data via a bend in the dispersion and an ac- 
companying change in spectral intensity, (b) Band dispersion 
derived by fitting the momentum distribution curve (MDCs, 
intensity as a function of momentum at fixed energy) at each 
energy in (a) to a Lorentzian. The low-energy kink is defined 
by the deviation of the dispersion from Vmid, the velocity fit 
between 30-40 meV (dotted fine), (c) MDC FWHM shows 
a more rapid decrease close to Ef as a consequence of the 
low-energy kink. 



of the high-Tc cuprates. A corresponding signature is 
also seen in the MDC FWHM, which shows a more rapid 
decrease for energies below the low-energy kink. Addi- 
tionally, these features have been observed over a wide 
doping range in underdoped Bi-2212, strengthening with 
underdoping, suggesting that it may be a ubiquitous as- 
pect of nodal physics.^ 

The 70 meV kink has been interpreted in terms of cou- 
pling to a sharp bosonic mode of energy J7 which is of 
either a lattice^^"^^ or magnetic^^"^^ origin. By anal- 
ogy, it has been suggested that the low-energy kink can 
be interpreted as coupling to an optical phonon.^ In Bi- 
2212 there are known infrared active optical modes at 97 
(12) and 117 cm~^ (14.5 meV) involving c-axis motion 
of the Cu,Ca,Sr and Bi ions^^ as well as a Raman ac- 
tive c-axis optical branch between 58 and 65 cm~^ (7.2 
and 8.1 meV) with the latter recently invoked to explain 
the low-energy kink.^ However, coupling to these optical 
modes is likely to be electrostatic in nature and therefore 
not sharply peaked in momentum space. Such a coupling 
would would produce a kink at $7 + Aq.^^"^^'^*^ However, 
it is possible to have a non-gap-shifted feature in the self- 
energy near the node if the coupling is strongly peaked 
in the forward scattering direction. Such is the case for 
the in-plane acoustic mode, which we consider here. 

The generic form for the el-ph coupling Hamiltonian is 



given by 

k,q,o',i.' 

(1) 

where dj^ g. (dk.o-) creates (annihilates) an electron in 
the antibonding pd — a* band with momentum k and 
spin tj and ^ ^<a,v) creates (annihilates) a phonon 
quanta of momentum q in branch v. For the in- 
plane acoustic branch, the el-ph interaction arises via a 
deformation-type coupling of the periodic lattice poten- 
tial. In this case, momentum-dependent el-ph coupling 
constant ^(k, q) is only a function of q and is given by^^ 

/ ^ 1 \ ft - ^(q) 

^(^^=V^\/2AflXq)^'>-^^ 

with fi(q) the phonon dispersion, e(q) = 1 + q^pjq^ 
the Thomas- Fermi dielectric function, q^p the Thomas- 
Fermi wavevector, Vceii = o,^c the unit cell volume, M the 
copper -I- oxygen ion mass, Cq the phonon polarization 
vector and V{q) — Aire'^ /q^e the Coulomb potential with 
e the static dielectric constant. 

Apart from the dispersion of il(q), Eq. (2) has a q- 
dependence governed by cx Cq • q/(q^ + Qtf) which is 
strongly peaked for q ~ qxF but with ^(q) — )■ for q 0. 
In the underdoped cuprates the Thomas-Fermi wavevec- 
tor is small and |ff(q)P is sharply peaked for a subset of 
small q.^'' As a result, the el-ph self-energy is dominated 
by contributions from scattering process from states i?k 
to nearby states E'k+q, with E^{k) = ^^(k) + A'^{k), 
where ^(k) = e(k) — /i is the electron dispersion in the 
normal state measured relative to the Fermi level. The 
nodal self-energy is thus determined from scattering to 
nearby states with a small superconducting gap and thus 
the peak in the self-energy is not shifted by the full gap 
but rather an average of the gap near the node. 

To demonstrate this we now calculate the spectral 
function for coupling to the acoustic phonon branch with 

0(q) = rJoysii?(g|,|a/2yTsin^^(jj/2y2)/^ and = 15 

meV.^^ We also include coupling to the 55 and 36 meV 
in- and out-of-phase c-axis polarized modes (the so-called 
Aig and Big branches, respectively), as well as the 70 
meV Cu-0 bond-stretching modes in order to also cap- 
ture the well-known 70 meV kink and renormalizations at 
higher binding energy. Since our focus is on the features 
at low energy, we treat the optical modes as dispersion- 
less with a momentum-independent coupling. The su- 
perconducting gap is modelled with a pure d-wavc form 
A(k) = Ao[cos(A:a;a) — cos{kya)]/2, where Aq = 37 meV 
is the maximum value at (0, tt). We take the dielectric 
constant to be set by a large in-plane value with e = 30eo 
(in this context we regard e as a free parameter set 
to obtain overall agreement with experiment) and take 
Qtf — O-bir/a. Finally, we take the effective lattice con- 
stants per plane to be a = 6 = 3.8 and c = 7.65 A, appro- 
priate for Bi-2212. The self-energy E(k, w) and spectral 
function A(k, oj) are calculated within Migdal-Eliashberg 
theory, which is the same formalism used in prior work to 
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FIG. 2: (Color online) (a) The calculated spectral function 
yl(k, Lj) along the nodal direction (0, 0)-(7r, tt) of supercon- 
ducting (25 K) Bi-2212. (b) The MDC-derived dispersion 
and (c) MDC FWHM obtained from the calculated A(k,cj). 
In calculating A(k,ij), an additional component cx ijp" has 
been added to the imaginary part of the self-energy in order 
to simulate additional broadening due to el-el interactions. 



understand the higher energy renormaUzations^'^^^'*'^"^. In 
order to introduce increased broadening in the spectral 
function with underdoping, we evaluate the self-energies 
while adding a finite imaginary part (b — 1 meV) to the 
bare electron Green's function. 

The strength of the el-ph coupling can be estimated 
from the slope of the self-energy at the Fermi level 



A(k) 



9ReI](k,w) 



(3) 



LJ=0 



For the above choice of parameters we obtain A(k„ode) = 
0.56 for the acoustic branch, consistent with the ob- 
served magnitude of the renormalization.^''^'^ For the opti- 
cal modes we set the overall strength such that Az — 0.23, 
0.05 and 0.05 for the 36, 55 and 70 meV modes, respec- 
tively, in order to match previous work.^'^'^^ 

The calculated ^(k, uj) (and corresponding MDC anal- 
ysis) for nodal cut in the superconducting state (25 K) is 
shown in Fig. 2. As with the experimental data (Fig. 1), 
the low energy kink is evident in the raw spectral func- 
tion appearing at an energy set by the phonon dispersion 
weighted by g(q). Since the coupling for scattering to 
the antinode is suppressed, this energy scale is slightly 
smaller than the van Hove points in the phonon density 
of states (10.6 and 15 meV). The low-energy kink is also 
evident in the MDC-derived dispersion and linewidth, 
shown in Figs. 2b and 2c, respectively. Both are in ex- 
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FIG. 3: Nodal MDC derived dispersions obtained from 
the model j4(k, lo) for various values of the Thomas-Fermi 
wavevector. The value of the dielectric constant e has been 
held fixed for each case. 



cellent agreement with the experimental data (Fig. 1), 
demonstrating that the forward scattering nature of the 
coupling to the acoustic branch accounts for both the en- 
ergy scale of the low-energy kink and the overall strength 
of the renormalization near the Fermi level. 

Since the coupling to the acoustic mode arises from 
the modulation of the screened Coulomb interaction a 
doping dependence naturally emerges from the increased 
metallicity of the cuprates with overdoping.^"^ To exam- 
ine the implications of this, in Fig. 3 we plot nodal 
MDC-derived dispersions obtained for our model calcu- 
lations while varying q^F from 0.47r/a to 0.67r/a. This 
allows us to mimic the increased screening which occurs 
with progressive overdoping. For decreasing qxF the 
overall strength of the coupling increases and becomes 
increasingly peaked at small q. As a result, stronger 
kinks are produced with underdoping in agreement with 
experiment.^ We also note that increased coupling at 
small q with underdoping has also been reported for 
more a complicated form for the dielectric function based 
on RPA treatments'^ and is therefore expected to be a 
generic result. These effects are also expected to sur- 
vive the inclusion of the short-range Hubbard interaction, 
which largely suppresses the el-ph vertex at large q over 
its value at small q.'^ 

We now turn our attention to the momentum depen- 
dence of the low-energy kink. As previously mentioned, 
for small qTP the el-ph self-energy is dominated by scat- 
tering to nearby momentum states and thus samples the 
local k-space value of the superconducting gap. There- 
fore, for off-nodal cuts, the energy of the kink should 
begin to gap shift, following the local value of A(k) for 
each cut. In Figs. 4a and 4b we present experimental 
and calculated MDC-derived dispersions for a series of 
off-nodal cuts in a direction parallel to (0,0)-(7r, tt) with 
cuts taken up to 11 degrees away from the node. Ex- 
perimentally, the superconducting gap is determined in- 
dependently for each cut by fitting symmetrized energy 
distribution curves at kp to a minimal model proposed 
by Norman et al.?'^ In Fig. 4a the kink position is de- 
termined from the deviation of the a straight line de- 
termined from fitting the dispersion slope v„iid between 
30-40 meV. The parameters for the calculations (Fig. 4b) 
are identical to those used in Fig. 2. 
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FIG. 4: (Color online) Low energy kink away from node (UD55, IIK). (a) MDC dispersions at node (rightmost) and away from 
the node, offset horizontally for clarity. Black circle represents gap energy, which is determined by fitting symmetrized data to a 
minimal model. Red circle defines the approximate low-energy kink position, as determined by where the measured dispersion 
deviates from Vmid- (b) MDC dispersions obtained from the calculated A(k, a;), (c) A sketch of how the kink position can be 
further quantified by subtracting a linear offset from the dispersion (shown here for nodal cut) between the gap energy and 
40 meV. (d) Difference in k between dispersion and dotted line in (c) for cuts from the node (top) and moving away from the 
node (towards the bottom). The data has is smoothed over 20 iterations and the curves have been offset vertically for clarity. 
Black dots indicate the gap energies and arrows mark peaks of these curves, determined from derivative after smoothing, (e) 
The difference in k between the MDC dispersion and linear dispersion as in (d) but obtained from the calculated j4(k,aj). (f) 
Momentum dependence of gap and kink energies, the latter determined by method in (c). 



As can be seen in the experimental data, as cuts are 
taken progressively further from the node both the su- 
perconducting gap and kink position shift to higher en- 
ergies. This behavior is also reflected in the linear bare 
band subtraction, shown in Figs. 4d and is qualitatively 
reproduced by the model calculations. The momentum 
dependence of both energy scales are summarized in Fig. 
4f, vifhich shows a clear, constant offset between the lo- 
cal value of the gap and the kink position, similar to 
a previous study. ^ Fig. 4b shows the calculated MDC 
dispersions and Fig. 4e show the subtracted linear bare 
band (Fig. 4c) obtained from applying the same proce- 
dures to the calculated spectral functions. The agree- 
ment between the theory and experiment is excellent 
and show that a strong coupling to the in-plane acoustic 
branch also captures the momentum-dependence of the 
low-energy kink. 

Our results indicate that the acoustic modes cou- 
ple to the lattice in a strongly doping and momentum- 
dependent way, which have a number of interesting im- 
plications and general consequences for these and other 
phonon modes. As the Coulomb interaction may dramat- 
ically alter the momentum-dependence of lattice effects, 
which become more pronounced as the screening length 
increases, the effect on transport can be substantial.^"^ 
Due to the large forward scattering peak, coupling to the 



acoustic branch may become less prominent in transport 
measurements apart from its overall renormalization of 
the Fermi velocity, as noted in a previous study (Ref. 9). 
Importantly, the d^2_y2 pairing interaction from phonons 
can be dramatically changed whereby the screened e-ph 
interaction which favors large momentum scattering can 
become more peaked at small q as screening breaks down 
and therefore might not as detrimental to d-wave pairing. 
This implies that the acoustic modes may in fact help d- 
wave pairing, in conjunction with electrons' coupling to 
c-axis and bond-stretching modes, in a way that would 
not be captured in fully itinerant and weak coupling sce- 
narios. This further highlights the need to consider lat- 
tice effects and coupling to the full spectrum of oxygen 
phonon modes in order to obtain a complete understand- 
ing of the the physics of the high-Tc cuprates. 
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